Goldstone boson counting in linear sigma models with chemical potential 



Tomas BrauneiQ 

Department of Theoretical Physics, Nuclear Physics Institute, 25068 Rez, Czech Republic 

We analyze the effects of finite chemical potential on spontaneous breaking of internal symmetries 
within the class of relativistic field theories described by the linear sigma model. Special attention is 
paid to the emergence of "abnormal" Goldstone bosons with quadratic dispersion relation. We show 
that their presence is tightly connected to nonzero density of the Noether charges, and formulate 
a general counting rule. The general results are demonstrated on an SU(3) x U(l) invariant model 
with an SU(3)-sextet scalar field, which describes one of the color-superconducting phases of QCD. 
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I. INTRODUCTION 

Spontaneous symmetry breaking plays an important 
role in many areas of physics and encounters a host of 
fascinating phenomena. The most distinguishing feature 
of spontaneous symmetry breaking is the presence of soft 
modes, long-wavelength fluctuations of the order param- 
eters), guaranteed by the Goldstone theorem [HQ. 

For low-energy properties of the spontaneously bro- 
ken symmetry it is important to know the number of 
the Goldstone bosons (GBs). While for spontaneously 
broken internal symmetry (space-time symmetries will 
not be the subject of this paper, see e.g. Ref. Q) 
in a Lorentz-invariant field theory it is always equal to 
the number of broken symmetry generators, the original 
Goldstone theorem predicts the existence of at least one 
GB. Indeed, there are several examples in nonrelativis- 
tic physics where the number of GBs is smaller than one 
would naively expect. The most profound one is perhaps 
the ferromagnet where the rotational SO (3) symmetry is 
spontaneously broken down to SO (2), but only one GB 
(the magnon) exists. 

The issue of GB counting in nonrelativistic field the- 
ories was enlightened by Nielsen and Chadha They 
showed that the defect in the number of GBs is related 
to the low-momentum behavior of their dispersion rela- 
tions. GBs with energy proportional to an odd power of 
momentum are classified as type-I, and those with energy 
proportional to an even power of momentum as type-II. 
The improved counting rule then states that the number 
of GBs of type I plus twice the number of GBs of type II 
is greater or equal to the number of broken generators. 

It should be noted that the form of the disper- 
sion law of the lightest degrees of freedom has impor- 
tant phenomenological consequences, e.g. for the low- 
temperature thermodynamics of the system. For in- 
stance, the heat capacity of a gas of bosons with E oc |p| 
falls down as T 3 for T — ► 0, while for bosons with i?ocp 2 
it is only T 3 / 2 . If no massless particles are present, the 
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heat capacity is suppressed by factor e m / kT : where m 
is the mass of the lightest particle. 

The interest in the problem of GB counting has been 
revived recently, mainly thanks to the progress in un- 
derstanding the phase diagram of quantum chromody- 
namics. At finite density Lorentz invariance is explicitly 
broken and GBs with nonlinear (as a matter of fact, gen- 
erally quadratic) dispersion relations may appear even in 
a relativistic field theory as a medium effect j^, [(j . Their 
presence turns out to be connected to the fact that some 
of the broken Noether charges develop nonzero density in 
the ground state, as has been observed in various color- 
superconducting phases of QCD [3, |8| or in a neutron 
ferromagnet 0- 

Schafer et al. have proved the following theorem: 
if the commutators of all pairs of broken generators have 
zero ground-state expectation value, then the number of 
GBs is equal to the number of broken generators. It 
is therefore clear that the nonzero charge density itself 
is not sufficient for a quadratic GB to appear. Indeed, 
the baryon number density does not make any harm to 
the usual linear GBs in the color superconductors. The 
corresponding generator must rather be a part of a non- 
Abelian symmetry group. Our main goal is to show that 
the opposite to the theorem of Schafer et al. generally 
holds: nonzero density of a commutator of two broken 
generators implies one GB with quadratic dispersion law. 

The paper is organized as follows. The following sec- 
tion is devoted to preparatory considerations: we explain 
how the quadratic GB is manifested in the Goldstone 
commutator and sketch its realization in the linear sigma 
model. In the next part, an example with an SU(3)- 
sextet condensation is investigated in detail. The general 
analysis is performed in the last section. 



II. PRELIMINARY CONSIDERATIONS 

In this section we shall investigate how the quadratic 
GBs come about, first at the rather general level of the 
Goldstone commutator and later more explicitly within 
the linear sigma model. 
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A. Goldstone commutator 

Let us briefly recall the proof of the Goldstone theorem. 
Following Ref. p|, we assume there is a local (possibly 
composite) field &(x) and a broken Noether charge Q 
such that (0|[<l>(x), Q]\0) ^ 0. Inserting the complete set 
of intermediate states into the commutator, one arrives 
at the representation 



<0|[4(aO,Q]|0) = [e^ Bkt (0|<i>(0)|n k )(n k |j o (0)|0) 



n=l 



iE-kt /n| „-0 



(0|j°(0)|n_ k )<n_ k |$(0)|0)] at k = 0, (1) 

where the index n counts the GBs. 

Now assume that we deal with a non-Abelian sym- 
metry group and some of its charges have nonzero den- 
sity in the ground state. Take as the GB field &(x) 
the zero component of the Noether current itself, so 
that {0\bi(x),Q b ]\Q) = */a4c(0U°(x)|0>, where f abc is 
the set of structure constants of the symmetry group. 
Should this be nonzero, we infer from Eq. Q that both 
(0|j°(0)|n) and (n\j® (0)|0) must be nonzero for some 
Goldstone mode n. 

The point of the above heuristic argument is that while 
in Lorentz-invariant theories there is a one-to-one cor- 
respondence between the GBs and the broken currents, 
here a single GB couples to two Noether currents. This 
explains (not proves, of course) at a very elementary level 
how the GB counting rule is to be modified in the pres- 
ence of nonzero charge density. 

One should perhaps note that the Nielsen-Chadha 
counting rule is formulated in terms of the GB disper- 
sion relations rather then charge densities. The connec- 
tion between these two was clarified by Leutwyler |10| . 
who showed by the analysis of the Ward identities for 
the broken symmetry, that nonzero density of a non- 
Abelian charge induces a term in the low-energy effective 
Lagrangian with a single time derivative. The leading or- 
der effective Lagrangian is thus of the Schrodinger type 
and the energy of the GB is proportional to momentum 
squared. 



B. Goldstone bosons within the linear sigma model 

In order to elaborate more on the properties of the 
GBs, we restrict ourselves from now on to the framework 
of the linear sigma model, that is a general scalar field 
theory with quartic self-interaction. 

To see how the Goldstone commutator emerges in this 
language, recall the SU(2) x U(l) invariant model of 
Schafer et al. [5j and Miransky and Shovkovy 6]. The 
Lagrangian for the complex doublet field of mass M in 
Minkowski space reads 



Finite density of the statistical system is represented by 
the chemical potential /i, which enters the Lagrangian 
in terms of the covariant derivative D^cj) — (dp — 
i$0ii^)4>- Upon expanding the covariant derivatives, the 
Lagrangian becomes 

C = d^Wf- 2nlm^do4>+ (^ 2 - M 2 )$(t>- A(<^) 2 . 

(2) 

For fj, > M the static potential develops a nontrivial 
minimum and the scalar field condenses. To find the 
spectrum of excitations at tree level we reparameterize it 
as 
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__ c in k T k /v 

V2 \v + ip 







H 2 -M 2 
A 



and look at the bilinear part of the Lagrangian. The 
crucial contribution comes from the term in Eq. J3J) with 
one time derivative. Upon expanding the exponentials it 
yields among others the expression 



--Him ( 1 ) [ir k T k , doirm] ( ^ 



= ^(7ri9o7r 2 -7r 2 9o7ri)- 



As will be made clear in the next subsection, it is this 
term that is responsible for the quadratic dispersion re- 
lation of one of the GBs. Its origin from the nonzero 
density of a commutator of two generators is now made 
obvious. This is the main idea to be remembered. The 
necessary technical details will come in the next two sec- 
tions. 



C. Bilinear Lagrangians and dispersion laws 

Bilinear Lagrangians with single-time-derivative terms 
will frequently occur throughout the whole text. It is 
therefore worthwhile to fix once for all the corresponding 
excitation spectrum. 

The bilinear Lagrangians we will encounter will have 
the generic form 



1 



\{d,Hf - \f 2 {n)H 2 



g(fi)Hdon. 

(3) 

The notation suggests that H is a massive (Higgs) mode 
whose mass function / 2 (/x) depends on the chemical po- 
tential, while 7r is the Goldstone mode. The excitation 
spectrum is found from the poles of the two-point Green 
functions or, equivalently, by solving the condition 



det 



E 2 p 2 

-iEg{p) E 2 



-iEg(ii) 



0. 



It turns out there is one massive mode, with dispersion 
relation 



^ 2 -/ 2 (m)+.9 2 (m) + 0(p 2 ), 



(4) 
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and one massless mode, with dispersion relation 



E 2 = 



/ 2 (M) 



p(v) + g 2 (v) [/ 2 (m) + 9 2 (m)] c 



o(p 6 ). 



(5) 

Now if / 2 (^) > 0, the Lagrangian J3J indeed describes 
a massive particle and a GB, whose energy is linear in 
momentum in the long- wavelength limit. On the other 
hand, when / 2 (/i) = 0, that is when both it and H would 
correspond to linear GBs in the absence of the chemi- 
cal potential, the dispersion relation of the gapless mode 
reduces to E = p 2 /|g(^)|. This is the sought quadratic 
Goldstone. 

In conclusion, the term with a single time derivative in 
general mixes the original fields in the Lagrangian. Mix- 
ing of a massive mode with a massless one yields one 
massive particle and one linear GB, mixing of two mass- 
less modes results in a massive particle and a quadratic 
GB jig. 



III. LINEAR SIGMA MODEL FOR 
SU(3)-SEXTET CONDENSATION 

As a nontrivial demonstration of the general idea pro- 
posed in the previous section, we shall now analyze 
in detail a particular model of spontaneous symmetry 
breaking. Consider a scalar field <E> that transforms as 
a symmetric rank-two tensor under the group SU(3), 
$ — > UQU T . Such a field describes a one- flavor diquark 
condensate in one of the superconducting phases of QCD 

m 

In addition to the SU(3) group, <I> is subject to U(l) 
transformations corresponding to quark number, $ — > 
e ie $e ie = e 2i6 <&. The most general SU(3)xU(l) invariant 
Lagrangian has the form 

C = tr(L> Al $ t L>' I $)-M 2 tr$ t <I>-atr($ t $) 2 -fe(tr$ t $) 2 . 

The quark- number U(l) has been assigned chemical po- 
tential fj, so that D & = (do — 2i^)$. The parameters 
a, b are constrained by the requirement of boundedness of 
the static potential [12{ . It is necessary that either both 
are non-negative (and at least one of them nonzero), or 
a < and b > \a\, or b < and a > 3\b\. 



A. Minimum of the static potential 

We start our analysis with a careful inspection of the 
static potential, 

= -(4/i 2 - M 2 ) tr $ f $ + atr($ t $) 2 + 6(tr $ f $) 2 . 

(7) 

A potential of the same type has been analyzed by Iida 
and Baym 0]. In their case, however, the global sym- 
metry was different, and we therefore provide full details. 

When A\i 2 - M 2 > 0, the stationary point <f> = be- 
comes unstable and a new, nontrivial minimum appears 



|2fj. The stationary-point condition reads 

$(-4/i 2 + Af 2 + 2a$ t $ + 2&tr$ t $) =0. (8) 

Before going into detailed solution of this equation we 
note that by multiplying Eq. JSJ from left by & and tak- 
ing the trace, the stationary-point value of the potential 
Q) is found to be 



-(4/i 2 - M 2 )tr$ f $. 



Any nontrivial stationary point of the potential is thus 
energetically more favorable than the perturbative vac- 
uum <i> = 0. We are, however, obliged to find a stable 
ground state, that is the absolute minimum of the poten- 
tial. 

We now make use of the fact that the field $ can always 
be brought by a suitable SU(3) x U(l) transformation to 
the standard form, which is a real diagonal matrix with 
non-negative entries ^4[. Eq. (|SJ then splits into three 
conditions and it is easy to see that all nonzero diagonal 
elements acquire the same value, denoted here by A. 

Let there be n of them, n = 1,2,3. Eq. implies 



A z = - 



1 A[x 2 - M 



2\2 



bn 



V s tnt — r 



1 {An 2 - M ) 



b+± 



To find the absolute minimum of the potential, it remains 
to minimize this expression with respect to n. 

For a > the minimum occurs at n = 3, and $ is 
proportional to the unit matrix, $ = Al, where 



1 4 M 2 
2 a- 



M 2 



36 



The SU(3) x U(l) symmetry is broken down to SO(3). 

For a < the potential is minimized by n = 1, that is 
$ is diagonal with a single nonzero entry and is conven- 
tionally chosen to be $ = diag(0, 0, A), where now 



M 2 



2 a 



The unbroken subgroup is now SU(2) x U(l). 

For a = the local minima corresponding to different 
n are degenerate since in that case, the Lagrangian © 
is invariant under an enhanced SU(6) x U(l) symmetry, 
treating $ as a fundamental sextet. Nonzero ground- 
state expectation value of <f> breaks this symmetry to 
SU(5) x U(l). As we shall see, such an enhanced symme- 
try leads to an increased number of GBs with quadratic 
dispersion relation [Til ]. 



B. Noether currents and charge densities 

Having found the vacuum configuration of the scalar 
field, we are ready to reparameterize it and find the exci- 
tation spectrum from the bilinear part of the Lagrangian. 
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Before doing that, we evaluate the ground-state densities 
of the Noether charges in order to make a priori predic- 
tions about the nature of the GBs. 

The infinitesimal SU(3) x U(l) transformation of <& 
has the generic form 5$ = i8k(\k& + ^A^), where the 
Afe stands for the Gell-Mann matrices (k — 1, . . . , 8) and 
the unit matrix (k = 0), respectively. The corresponding 
Noether currents are 



Jk 



■ tr [Z^$t ( Afe $ + $xl) - h.c] . 



Taking a generic static field configuration to be <I> = 
diag(Ai, A 2 , A 3 ) results in the charge densities 



Jo 

J3 



MA 2 



A| 
A?) 



A|), 



is° = ^(A? 



A^-2A 2 ). 



In the SO (3) symmetric phase (a > 0), all generators 
but the U(l) quark number have zero density. As this 
is an Abelian generator, we expect six linear GBs corre- 
sponding to the six broken generators 1, Ai, A3, A4, Xq, Xg. 
In the a < case, the densities of Ao and As are nonzero. 
This means that the commutators [A^As] and [Ag,A7] 
have nonzero ground-state density. With regard to the 
general discussion above, we thus expect two quadratic 
GBs corresponding to pairs (A4, A5) and (Xq, A7), and one 
linear GB of the generator As- 



C. The a > case 

We shall now proceed to the calculation of the mass 
spectrum of the a > phase. We could do well with 
just shifting $ by its vacuum expectation value Al, but 
this would complicate the identification of the massless 
modes. It is more convenient, and physical, to find such a 
parameterization that the GBs disappear from the static 
potential. 

To that end, recall that the field (now coordinate- 
dependent) can be brought to the diagonal form by a 
suitable SU(3) x U(l) transformation. In other words, it 
may be written as 

$(a;) = e 2l6 ^U(x)D(x)U T (x), 

where U(x) £ SU(3) and D(x) is real, diagonal and 
non-negative. Now the unitary matrix U can be (at 
least in the vicinity of unity) expressed as a product 
U = VO, O € SO(3) being an element of the unbroken 
subgroup and V being built from the broken generators, 
V = e^k^k^ k = 1,3,4,6,8. A simple observation that 
0(x)D(x)0 T (x) is the general parameterization of a real 
symmetric matrix leads to the final prescription, 

= e 2ie ^V(x) [Al + ip(x)} V T (x). 

The real symmetric matrix ip contains six massive modes, 
while V contains five GBs. With 8 this is altogether 



twelve degrees of freedom, as it should for $ is a complex 
symmetric 3x3 matrix. 

It is now straightforward, though somewhat tedious, to 
plug this parameterization into the Lagrangian JSJ and 
expand to the second order in the fields. Omitting details 
of the calculations, we just report on the results. 

The full static potential (up to a constant term - the 
vacuum energy density) becomes 

V($) = 4A 2 [a tr tp 2 + 6(tr Lpf] 

+ 4A (atrip 3 + btriptrip 2 ) + atrtp 4 + b(tnp 2 ) 2 . 

The bilinear Lagrangian turns out to be (we use the no- 
tation V = e ,;n ) 

£ bili „ = 12 A 2 {d^ef + 4A 2 tr(^n) 2 + tr(d^) 2 
-4A 2 [a tr ip 2 + 6(tr ip) 2 ] -16^A [d 6 tr ip + tr(^<9 II)] . 

The kinetic terms are brought to the canonical form by 
a simple rescaling of the fields, upon which the spectrum 
is readily determined from Eqs. I@J and ©. 

The excitations fall into irreducible multiplets of the 
unbroken SO(3) group. There are two singlets, stemming 
from the mixing of 9 and tr ip, 

massive mode E 2 = 24/i 2 - 2Af 2 + C(p 2 ), 
linear GB E 2 = ^ p 2 + 0( P % 

and two 5-plets, the mixtures of 713, 7r4, ttq, tts) and 
the traceless part of tp, 



massive modes E 



a _ (24 M 2 -2M> + 48/i 2 b + 2 
a + 36 

linear GBs E 2 = —pl—M\—p 2 + 0(p 4 ). 

(12/i 2 - M 2 )a + 24/i 2 6 v ' 

It is easily seen from these formulas that the masses of 
the massive singlet and the massive 5-plet are connected 
by 

ml =ml + (An 2 - M 2 ) = ml + 12A 2 6. 

a + 3b 

The singlet is heavier than the 5-plet for b > and vice 
versa. 

The excitation spectrum is plotted in Fig. ^ for the 
case M 2 > 0. Below the phase transition to the Bose- 
Einstein-condensed phase, the medium-modified disper- 
sion relations are simply E — \/p 2 + M 2 ± 2/i. Right at 
the transition point, there are six modes with mass 2M 
and six massless ones with dispersion E = p 2 /4:(i. As the 
phase transition is second order, the dispersion relations 
of all excitation branches must be continuous functions 
of /i, that is all GBs become quadratic at the transition 
point. This is also easily checked on the broken-symmetry 
side of the transition. As 2fi — > M+, the phase veloci- 
ties of the linear GBs tend to zero, and their dispersions 
become quadratic. 
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FIG. 1: Mass spectrum as a function of the chemical potential 
for a > 0. The boldface-typed numbers denote the degenera- 
cies of the excitation branches. To obtain numerical results, 
particular values a = b = 1 were chosen. 



Note that also for a = the dispersion relation of the 
GB 5-plet becomes quadratic, E — p 2 /4fi. This is in 
accord with the enhanced SU(6) x U(l) symmetry of the 
Lagrangian. There are altogether eleven broken gener- 
ators of the coset SU(6)/SU(5), one linear GB and five 
quadratic ones [forming now the 5-plet of the unbroken 
SU(5)], and the Nielsen-Chadha counting rule is thus 
satisfied. 



D. The a < case 

We use the same method for parameterization of <& 
as in the previous case. This time we write $(x) = 
U(x)D(x)U T (x), where U{x) £ SU(3) x U(l). Next per- 
form the decomposition U = e lII U', where n = irk^k, 
k = 4, 5, 6, 7, 8, and U' belongs to the unbroken subgroup 
SU(2) x U(l). Since U' (x)D(x)U' T (x) is block-diagonal 
with a complex symmetric 2x2 matrix in the upper-left 
corner, we arrive at the parameterization 



= e m{x) [diag(0, 0, A) + E(a;)] e mT(x \ E(x) = 









H{x) t 



Here H is a real field and a is a complex symmetric 2x2 matrix. These two embody the massive modes that survive 
in the static potential, 

V ($) = 4A 2 (a + b)H 2 - 2A 2 a tr ctV + 4A(a + b)H 3 + (a + b)H 4 + AAbH tr a* a + 2bH 2 tr <Ao + a tr(a^a) 2 + 6(tr crV) 2 . 
The bilinear part of the Lagrangian reads 

-Cbilin = tr(d^d^a) + {d^H) 2 + 2A 2 (9 M n^n) 33 + 2A 2 (a M n 33 ) 2 - 4A 2 (a + b)H 2 + 2A 2 atr ( r t a 

- 16/iA7?a n 33 - 4/xA 2 Im[n,don] 33 - 4/zImtrcr t S cr. 



The excitations are again organized in multiplets of 
the unbroken SU(2) x U(l). H and 7r 8 mix to form two 
singlets, 

massive mode E 2 = 24^i 2 - 2Af 2 + C(p 2 ), 
linear GB E* = ^^ P 2 + 0(p% 

and the pairs (774, tts) and (tt6, ^7) give rise to a doublet 
of massive modes and a doublet of massless ones, 

massive modes E 2 = 16^ 2 + C(p 2 ), 

D 4 

quadratic GBs E 2 = + Ofp 6 ). 

The matrix a represents two triplets of massive par- 
ticles. The part of the bilinear Lagrangian containing a 



may be rewritten as 

C a = triD^D^a) - (4^ 2 + 2A 2 |a|) trtrV, 
which immediately implies the dispersion relations 
E = v/4^ 2 + 2A 2 |a| ± 2fi + 0(p 2 ). 

The mass spectrum is shown in Fig. El The unbroken- 
phase part of the spectrum is the same as in the a > 
case, since for 2/x < M the tree-level masses of the 
particles do not depend at all on the quartic potential, 
i.e. the parameters a, b. Also, the same remark about 
the continuity of the dispersion relations across the phase 
transition applies. 

Again, in the limit a = 0, the lighter of the two triplets 
in a becomes a triplet of quadratic GBs, and joins the 
other two quadratic GBs to form the full SU(5) 5-plet. 
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FIG. 2: Mass spectrum as a function of the chemical potential 
for a < 0. The singlet and triplet lines are so close that they 
almost coincide, but they are not degenerate. The spectrum 
is plotted for a = —0.5 and 6=1. 



6 

5 



2 











\ . 


1 


1 




\ 2 


5 


V 




3^ 


6 



-1 



a 



FIG. 4: Mass spectrum as a function of a. The graph is 
plotted for fj, = M and 6 = 1. The potential is unstable for 
a < — 1. The singular behavior of the masses of a is due to 
divergence of A towards the stability limit of the potential. 




FIG. 3: Phase diagram of the linear sigma model for SU(3)- 
sextet condensation. The phases are labeled by the symmetry 
of the ground state. The line of second order phase transition 
at a = 0, b > has SU(5) x U(l) symmetry. 



To summarize our results, the theory described by the 
Lagrangian JBJ) has two different ordered phases, both 
occurring at 4/l/ 2 > M 2 , distinguished by the symmetry 
of the ground state. The corresponding phase diagram 
in the (a, b) plane is displayed in Fig. |21 

As the excitations above the ordered ground state are 
grouped into irreducible multiplets of the unbroken sym- 
metry, it is interesting to find out how the structure of 
these multiplets changes across the phase transition from 
one ordered phase to the other. In Fig. Q]we show the 
dependence of the masses on the parameter a at constant 
chemical potential. The masses are continuous functions 
of a as the transition is second order. 



As a final remark we note that in the original appli- 
cation of Ref. the field <& represented a diquark 
condensate and the SU(3) was the color gauge group of 
QCD. One might wonder whether the usual Higgs mech- 
anism for gauge boson masses survives when there are 
fewer GBs than the number of broken generators, be- 
cause of the presence of quadratic GBs. This question 
was answered affirmatively by Gusynin et al. jig, and 
there is therefore no need to worry about the fate of glu- 
ons. 



IV. GENERAL ANALYSIS 

In this section we shall collect experience gained by 
solving particular examples and set out for a general anal- 
ysis. We will find out, with some effort, that the ideas 
sketched in Sec. II and demonstrated in Sec. Ill have a 
straightforward generalization to a whole class of theo- 
ries. It is understood, however, that we shall all the time 
stay in the framework of the linear sigma model, and at 
the tree level. The possibilities of further progress are 
discussed in the conclusions. 



A. Chemical potential and global symmetry 

As the starting point we shall address the question 
what is the most general symmetry of a theory with 
nonzero chemical potential. 

Let the microscopic theory possess a global continu- 
ous symmetry with the corresponding conserved Noethcr 
charges. The physical meaning of the chemical potential 
\x is that we wish to fix the statistical average of a con- 
served charge, say Q. This is technically achieved by 
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introducing the grandcanonical ensemble and replacing 
the microscopic Hamiltonian H with H — /iQ. 

It is now clear that by adding the chemical potential, 
we break explicitly all Noether charges that do not com- 
mute with Q. This is the technical realization of the 
physically intuitive fact that we cannot keep simultane- 
ously fixed the values of two noncommuting operators 
(i.e. incompatible observables) . 

This simple observation implies that, as far as exact 
symmetry is concerned, chemical potential is always as- 
signed to a generator that commutes with all others, that 
is to a U(l) factor of the exact global symmetry group. 

Of course, when the symmetry of the microscopic the- 
ory is non-Abelian, then adding of the chemical potential 
generally produces a number of approximately conserved 
charges (at least for small f\i) that generate approximate 
symmetries. These may also be spontaneously broken, 
resulting in the corresponding set of pseudo-Goldstone 
bosons. Throughout this paper we are, however, con- 
cerned only with true GBs, and therefore only the exact 
global symmetry will be considered. 

It is also interesting to find out how the Abelian nature 
of the charge equipped with chemical potential is man- 
ifested in the Lagrangian formalism. There, as already 
mentioned, chemical potential enters the Lagrangian in 
terms of the covariant derivative of "matter" fields [Tl| . 

The Lagrangian can be made formally gauge-invariant 
by introducing an external gauge field A^. Provided 
the matter fields cf> transform under the symmetry group 
linearly as cf> — + U<f>, A^ transforms as usual as A^ — > 
UA^U^ 1 + iUdfJJ^ 1 . Now the exact symmetry is such 
that the Lagrangian is invariant under the global trans- 
formation of the matter fields with A^ fixed at A^ — 
(fxQ, 0,0,0), This is possible only when A^ = UA^U^ 1 . 
We thus again arrive at the conclusion that the generator 
being assigned chemical potential must commute with all 
others. 



B. Linear sigma model 

Now consider a general linear sigma model defined by 
the Lagrangian 

C = D^D»ct>-V{(t>). (9) 

Here <fi denotes a set of complex [2lJ scalar fields that 
form a (possibly reducible) multiplet of the exact global 
symmetry group G, i.e. span the target space of a (pos- 
sibly reducible) representation of G, say 1Z. V((f>) is 
the most general G-invariant static potential containing 
terms up to the fourth power of fa and the covariant 
derivative is given by D^cj) = (d^ —iA^ip. A^ is the con- 
stant external field that incorporates chemical potential 
for one or more U(l) factors of G, and is eventually set 
to Ap = (£\ faQm, 0, 0, 0), where the Qi's are the U(l) 
generators, the subscript 1Z denoting the image in the 
representation 1Z. 



Upon expanding the covariant derivatives Eq. @ 
takes the form 

C = d^d^ - 2 Im tfAfdrf - V eS (fa, (10) 

the effective /i-dependent potential being V e s(fa — 

Spontaneous symmetry breaking occurs when V c ff((f>) 
develops a nontrivial minimum at some <p = fa. In order 
to elucidate the physical content of such a theory, it is 
necessary to conveniently parameterize the field <f>. 

We stress the generality of the parameterization 
method suggested and applied in Sec. III. One first 
writes fax) = Uiz(x)fatd(x), where s td is a standard 
form to which the field <p can always be brought by a 
suitable transformation U € G. Next Uu is factorized as 
Uti = e m U^, where II is a linear combination of the bro- 
ken generators (or more precisely, their 72.-images) and U' 
belongs to the unbroken subgroup H. The final step is 
to identify UL (x)fatd(x) with a certain representation of 
H and parameterize it linearly as fa + H(x). H(x) is 
going to be the multiplet of massive (Higgs) fields. We 
therefore invoke the parameterization 

fax) = e m W [fa + H (x)} . (11) 

In order to specify the transformation properties of 
H, recall that the GBs transform linearly in the ad- 
joint representation of the unbroken subgroup |17| . i.e. 
n — > U^HU!^ 1 for any U' 6 H. As a consequence, 
H = e~ lIl (f) — fa transforms as H — * U^H, since fa 
is an H-singlet. 

To summarize, H transforms in the representation 1Z 
truncated to the subgroup H, and the multiplets of the 
massive modes are therefore found in the decomposition 
of 1Z into irreducible representations of H. 

For instance, in our case a > the symmetric rank- 
two tensor representation of SU(3) splits under the SO(3) 
subgroup into a traceless symmetric rank-two tensor and 
a singlet. On the other hand, in the a < case it yields 
a symmetric rank-two tensor of SU(2) (the field a) plus 
a singlet. 

As an aside let us remark that the physical spectrum 
of the theory of course does not depend on the param- 
eterization chosen for the field cf>. What if we chose e.g. 
the linear parameterization mentioned (and abandoned) 
above in Sec III C? Instead of Eq. (|llfl . we would then 
have analogously 

fax) = fa + H(x)+m(x)fa. (12) 

It is easy to see that the bilinear terms in the La- 
grangian with one or two derivatives come out identical 
as for the parameterization (|llf> . The reason is that the 
only difference stemming from the nonlinear structure of 
e lU could possibly come in the form fa\A^dfJI 2 fa, but 
this is real and therefore it drops out of the Lagrangian 

The only difficulty with the linear parameterization 
(|12fl is that the GBs do not disappear automatically from 
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the static potential. Instead, we have to use explicitly the 
G-invariance to show that II disappears from the bilinear 
(mass) part of the potential. 

Upon the field redefinition as in Eq. (|1 1|) , the effective 
potential V e g becomes (up to a constant term) 

V en (<f>) = V(<j)o + H) - (H^A^A^H + 2ReH^A^A^ ). 

As we are expanding the potential about its absolute min- 
imum, the additional term linear in H is right enough to 
cancel a similar term coming from V(4>o + H). We are 
interested in the bilinear part of the potential, Vhi\m{H) 1 
which determines the mass term for H. 

Now we analyze the first two terms of the Lagrangian 
(|IUII . The two-derivative term yields the bilinear contri- 
bution 

d^d^H + ^JfylWItyo + 2 Im (fad ^Jld^H. (13) 

The first two terms in Eq. (|I3[) are the expected kinetic 
terms for the Higgs and Goldstone fields, respectively. 
The GB term, however, asks for a check that it is non- 
degenerate. 

Let n(cc) = TTk(x)Tk, Tk being the set of broken genera- 
tors. The GB kinetic term becomes d^nkd^wi<fi TkTi(f>o = 
Id^TTkd^TTi^llTk, Ti}<?!> . The matrix 4>l{T k , T;}0 O is real 
and symmetric and may be chosen, by taking an appro- 
priate basis of broken generators, diagonal. It is obvi- 
ously nondegenerate, as necessary in order to have kinetic 
terms for all the GBs, since otherwise (fy^T^Tkipo — for 
some Tfc, implying that Tfc is in fact not broken. 

The third term in Eq. I|13[) eventually turns out to be 
zero. Nevertheless, as other terms of a similar structure 
will be dealt with in the following, we shall analyze it in 
detail. The crucial point is the way various fields trans- 
form under the unbroken subgroup H. Virtually all in- 
formation about the structure of the bilinear Lagrangian 
may be obtained by a proper decomposition of the rep- 
resentation 1Z into irreducible representations of H, and 
making repeated use of the Wigner-Eckart theorem. 

Now when H and II belong to different representa- 
tions of H, the Wigner-Eckart theorem immediately tells 
us that the last term of Eq. (|I3|) vanishes. There is, 
however, a subtle exception to this argument. As 1Z is 
a complex representation, real representations of H are 
doubled in its decomposition. The reason is that when 
the set of vectors Xfc constitute the basis of a real repre- 
sentation of H, the vectors i\k form an independent basis 
of an equivalent representation. 

It may be that H and n (or n^o) are such doubles. 
This happens, for instance, for the two 5-plets in Sec. Ill 
C. In such a case, however, (f^d^Hd^H is real and, again, 
does not contribute to Eq. (|I3|I . 

The single-derivative term in Eq. I|1(J|) gives, after a 
short manipulation, the bilinear terms 

-2 Im H^A^df.H-A Re H U^II0o-Im 0^ [II, d^IL^o. 

(14) 



Throughout the calculation we made use of the fact that 
A^ is a U(I) generator, and therefore commutes with n. 

Putting together all the pieces of Eqs. I|I3(I and l|I4|l . 
we arrive at our main result - the bilinear Lagrangian for 
a general linear sigma model, 

Cmm = d^d^H - VuxUH) - 2ImH^A^d ft H 

+4^na' i n(/>o-4 Re A^n^o-im 4>Ia^ [n, c^n]^. 

(15) 

This formula contains all information about the particle 
spectrum of the theory, and the rest of the section is 
therefore devoted to its analysis. 

C. Discussion of the results 

There are altogether three terms with a single time 
derivative in Eq. Qf 5p. The term Im_fft A^d^H causes 
splitting of the masses of the massive modes. The term 
Re A^d^Hipo mixes massive and massless modes and, 
according to Sec. II C produces linear GBs. Finally, the 
term Imc^^p, d^Ii]^ mixes the Goldstone fields and 
gives rise to the quadratic Goldstones. 

With the Wigner-Eckart theorem at hand it is easy 
to check that each of the elementary fields appears in at 
most one of the three single-derivative terms. This fact 
essentially reduces the analysis of the Lagrangian l|15|) to 
the model two-field problem discussed in Sec. II C. 

To prove it note that the mixing term Ke W A^d^Hipo 
can be nonzero only when H and n are the two copies 
of the doubled real representation of H. Now the real 
multiplet H gives real A^d^H , and therefore does 
not contribute to the mixing of the massive modes. 
The real multiplet n analogously does not contribute to 
Im^J^ITl, (9 M n]0o as a consequence of the analysis that 
follows. 

As the main concern of this paper is Goldstone boson 
counting, we shall now concentrate on the last term of 
Eq. (|I5|I . which produces the quadratic GBs. 

First, it is clear that our suspicion about the connec- 
tion between the quadratic GBs and nonzero charge den- 
sities was right. For by the very same method as in Sec. 
Ill B we derive the Noether current corresponding to the 
conserved charge T, 

jJ = -»(^-k), 

and the ground-state density of T is 

f T = 2^A 1> . (16) 

The last term of Eq. (|I5|I is therefore indeed proportional 
to the ground-state density of the commutator of two 
generators. 

We may now in the general case proceed as in Sec. 
Ill, that is find the ground state, calculate the Noether 
charge densities, and make a definite prediction for the 
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particle spectrum. We can, however, do even better, at 
least a bit. 

We need not calculate the charge densities explicitly 
to say, which of the generators may produce quadratic 
GBs. It is obvious from Eq. I|16fl that only such a gener- 
ator T may acquire nonzero density, which is a singlet of 
the unbroken subgroup H. We therefore just have to de- 
compose the adjoint representation of G into irreducible 
representations of H and look for spontaneously broken 
singlets. 

As in the examples above, we next choose such a ba- 
sis that all the generators with nonzero density mutually 
commute. This ensures that they can be completed to 
form the Cartan subalgebra of the Lie algebra of G. Fol- 
lowing the standard root decomposition of Lie algebras 
(see e.g. Ref. Qjj), the rest of generators group into pairs 
whose commutator lies in the Cartan subalgebra. They 
arc the lowering and raising operators or their hermitian 
linear combinations, and together with their commutator 
span an SU(2) subalgebra of G. 

The point of this procedure is that only pairs of Gold- 
stone fields are then mixed by the single-derivative term 
Im(/>jA M [n, (9 M n]</>o and the excitation spectrum may be 
fully described with the help of the simple two-field bi- 
linear Lagrangian 0. Consequently, the quadratic GBs 
count as one per each pair of generators whose commu- 
tator develops nonzero ground-state density. 

The feasibility of such a pairing also follows from group 
theory and the Wigner-Eckart theorem. As the commu- 
tator of the two generators is to be an H-singlet, they 
must come from the same irreducible representation of 
H. 

To briefly conclude this section, we once again empha- 
size the fact that almost all we need to know about the 
excitation spectrum of the general linear sigma model © 
may be extracted from the bilinear Lagrangian (|15|l by 
simple group theory. We decompose the adjoint repre- 
sentation of G with respect to the unbroken subgroup H 
to determine the multiplet structure of the Goldstones. 
The remaining H-multiplets in the decomposition of the 
representation 1Z of the scalar field <p are the massive 
modes. 

The quadratic GBs are discovered with the knowledge 
of the ground-state densities of the broken generators. 
Without further calculation, we can even determine their 
dispersion relations. Making use of the continuity of the 
dispersion relations across the phase transition and the 
known dispersion relations in the unbroken phase, we 
may assert that the quadratic GB dispersion relation is 
generically of the form E = p 2 /2/iQ, where Q is the 
charge of the GB field under the U(l) subgroup equipped 
with the chemical potential. 

V. CONCLUSIONS 

We have analyzed spontaneous breaking of internal 
symmetries in the framework of the relativistic linear 



sigma model with finite chemical potential. Our prime 
motivation was to establish a counting rule for Gold- 
stone bosons in view of the fact that explicit breaking 
of Lorentz invariance by medium effects may cause the 
number of GBs to differ from the number of broken sym- 
metry generators. 

Our results confirm the Nielsen-Chadha counting rule. 
We show that the GBs have either linear or quadratic 
dispersion law at low momentum, and that the number 
of the first plus twice the number of the second gives 
exactly the number of broken generators. 

In addition, we find a criterion which gives in a purely 
algebraic way the number of quadratic GBs, the only 
necessary input being the structure of the ground state. 
There is one quadratic GB for each pair of generators, 
whose commutator has nonzero ground-state density. 

However, despite the generality of our results, many 
open questions still remain. First, we stress the fact that 
we work all the time at the tree level. It would be inter- 
esting to know the effect of radiative corrections on the 
details of the spectrum. On the other hand, it seems that 
at least the dispersion relations of the quadratic GBs are 
rather generic as they depend only on the chemical poten- 
tial in a very simple way. There might be a more robust, 
nonperturbative method to determine them, which relies 
only on the broken symmetry, and does not depend on 
the details of the dynamics of symmetry breaking. 

Second, we worked within the linear sigma model as it 
is easy to manipulate perturbatively once the scalar field 
has been properly shifted to its new ground state. It may 
happen that our results are valid generally for relativistic 
theories with chemical potential. At least the argument 
presented in Sec. II A that clarifies the connection be- 
tween the charge densities and the GB counting, suggests 
such a possibility. 

As adding chemical potential breaks Lorentz invari- 
ance in a very particular way, it might be possible to 
strengthen the Nielsen-Chadha counting rule at the cost 
of limiting its validity to a smaller class of theories. Even 
such a theorem would, however, find many applications 
on relativistic many-particle systems. We hope that our 
future work will help to find the answer to these ques- 
tions. 
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